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Abstract 

Forward jet cross sections have been measured in neutral current deep inelastic 
scattering at low Bjorken-x with the ZEUS detector at HERA using an integrated 
luminosity of 81.8 pb~^. Measurements are presented for inclusive forward jets 
as well as for forward jets accompanied by a dijet system. The explored phase 
space, with jet pseudorapidity up to 4.3 is expected to be particularly sensitive 
to the dynamics of QCD parton evolution at low x. The measurements are 
compared to fixed-order QCD calculations and to leading-order parton-shower 
Monte Carlo models. 
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1 Introduction 



Deep inelastic lepton scattering (DIS) off protons provides a rich field for exploring the 
parton dynamics in QCD. HERA has extended the phase-space region in the Bjorken 
scaling variable, xbj, down to a few 10~^. At such low xbj, several steps in the QCD 
cascade initiated by a parton from the proton can occur before the final interaction with 
the virtual photon takes place. The result of this cascade may be observed in the final 
state and provides an opportunity to study the QCD parton evolution in detail. 

Within perturbative QCD (pQCD), fixed-order calculations for the parton evolution are 
so far available only at next-to- leading order (NLO). A number of different approximations 
to the QCD evolution have been developed, based on summing of particular subsets of 
diagrams in accordance with their importance in the phase space considered. 

The conventional DGLAP [1-3] approach sums up the leading logarithms in the virtuality 
of the exchanged boson, Q^, and is expected to be valid at not too small xbj and . At 
small XBj, a better approximation is expected to be provided by the BFKL formalism [4] 
which resums the leading logarithm terms in where x is the fractional longitudinal 
momentum of a parton. The CCFM [5,6] approach interpolates between the two types 
of evolution, DGLAP and BFKL. 

The DGLAP evolution equations have been successfully tested at HERA in inclusive 
measurements at low xbj and no indication of BFKL dynamics was observed. The dy- 
namics at low XBj can be further probed by measurements of the partonic final state 
that highlight the differences between predictions of the BFKL and DGLAP formalisms. 
BFKL evolution results in a larger fraction of small xbj events with forward jet^ than 
the DGLAP evolution. A forward jet is characterised by a high fractional longitudinal 
momentum, a;-''^* = where p is the proton momentum and 'p'^ is the longitudinal 

jet momentum [7]. 

A comparison of data on forward jets with the DGLAP leading-order parton-shower Monte 
Carlo programmes performed previously in DIS at HERA [8-10] has revealed a clear 
deficit of forward jets in the Monte Carlo. However, the addition of a parton cascade 
evolved according to DGLAP on the photon side has significantly improved the description 
[11]. The simulation based on the Color Dipole Model (CDM) [12-14], which includes 
parton emissions not ordered in transverse momentum, also succeeded in describing the 
data. Fixed-order NLO QCD calculations [15,16] were also compared to the forward jets 
measurements in more recent publications and failed to describe the data. These studies 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 
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were performed up to pseudorapidities 3 [17] and 2.7 [18]. 

In this paper, measurements of inclusive forward-jet cross sections for pseudorapidities of 
up to 4.3 are presented, based on a data sample which corresponds to a twofold increase 
in luminosity with respect to the previous ZEUS analysis [17]. Furthermore, a comparison 
of the measured cross sections with the Cascade Monte Carlo [19], based on the CCFM 
evolution, is presented. In addition, measurements of "forward jet + dijet" cross sections, 
as investigated by the HI Collaboration [18], are reported. These measurements explore 
parton evolution in a more exclusive way and are more sensitive to its details. 

2 Theoretical framework and phase-space definitions 

The DGLAP evolution equations, based on coUinear factorisation, assume that the dom- 
inant contribution to parton evolution comes from subsequent parton emissions that are 
strongly ordered in transverse momenta, fcr, where the parton with the largest in- 
teracts with the photon. In this formalism, only the leading terms in InQ^ in the QCD 
perturbative expansion are summed up. Since this approximation does not resum leading 
In 1/x terms, it may not be adequate at low Xbj • 

Contrary to the DGLAP approach, the BFKL evolution equation resums the leading 
Inl/x and the evolution proceeds over x at fixed Q^. The BFKL approach imposes no 
ordering in but strong ordering in x, with the low-x parton interacting with the photon. 
At small ^Bj, the difference between these approaches to the QCD parton evolution is 
expected to be most prominent for hard partons created at the beginning of the cascade, 
i.e. at pseudorapidities close to the proton (forward region). 

The CCFM approach is based on the idea of coherent gluon radiation, which leads to 
angular ordering of gluon emissions in the gluon ladder. It interpolates between the 
above two types of evolution, so it should be applicable over a larger phase-space region. 

A phenomenological approach to parton evolution is provided by CDM [12-14]. In this 
model, gluons are emitted by color dipoles successively spanned between partons in the 
cascade. Due to independent radiation of the dipoles, the emitted gluons are not ordered 
in kT and thus CDM mimics the BFKL-type evolution. 

To investigate the relevance of different approaches, events with a jet in the forward 
region were analysed in the Iow-xbj region. Events were required to have at least one jet 
satisfying the following criteria: 
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where p^^* is the transverse momentum of the jet. The first condition suppresses strong 
ordering in the transverse momenta and decreases the probabihty of having a DGLAP- 
type evolution. The second condition enhances the phase space for the BFKL evolution. 

A further event sample called "forward jet+dijet" , which contains at least two hard jets 
in addition to the forward jet (fjet), was selected. The jets were ordered in pseudorapidity 
such that rf'^^^ < if'^^^ < r]^^'^^. For this sample, the pseudorapidity separation of dijets, 
A?7i = T^j'^*^ — rf'^^^ , and the pseudorapidity difference between the forward and the second 
jet of the dijet, Ari2 = r/^-^*^* — 77-^°*^, were studied. 

The cross section as a function of A?72 was investigated for two intervals of Arji, Arji < 1 
and A?7i > 1. With such a choice of Arji, different dynamics of the partons in the cascade 
are expected to be highlighted. For Arii < 1, small invariant masses of the dijet system 
are favoured and, therefore, partons with small values of Xg are produced, where Xg is the 
longitudinal momentum fraction carried by the gluon coupled to the hard dijet system 
(Fig. [1]). Consequently, a large space is left for BFKL-type evolution in x from the forward 
jet to the dijet system. When Arji is large, BFKL-like evolution can occur between the 
partons producing the dijet system. 

3 Experimental set-up 

The analysis was performed with the data taken with the ZEUS detector from 1998 
to 2000, when HERA collided electrons or positron^ with energy of = 27.5 GeV with 
protons of energy Ep = 920 GeV, yielding a centre-of-mass energy of 318 GeV. The results 
are based on the sum of the e~p and e~^p samples, corresponding to integrated luminosities 
of 16.4 ± 0.3 pb~^ and 65.3 ± 1.5 pb~^, respectively. 

A detailed description of the ZEUS detector can be found elsewhere [20]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [21], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
polar-angle region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks is a{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [22] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 

^ Hereafter, both e+ and arc referred to as electrons, unless explicitly stated otherwise. 
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(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18/ \/E for electrons 
and a{E)/E = 0.35/ VE for hadrons, with E in GeV. 

The luminosity was measured using the bremsstrahlung process ep — e'-fp with the 
luminosity monitor [23], a lead-scintillator calorimeter placed in the HERA tunnel at 
Z = -107 m. 

For the 1998-2000 running period, the forward plug calorimeter (FPC) [24] was installed 
in the 20 x 20 cm^ beam hole of the FCAL, with a small hole of radius 3.15 cm in the center 
to accommodate the beam pipe. The FPC increased the forward calorimetric coverage 
by about 1 unit of pseudorapidity to t] <5. The FPC consisted of a lead-scintillator 
sandwich calorimeter divided longitudinally into electromagnetic and hadronic sections 
that were read out separately by wavelength-shifting fibers and photomultipliers. The 
energy resolution, as measured under test-beam conditions, was a{E)/E = 0.A1/\/E © 
0.062 and a(E)/E = 0.65/v^©0.06 for electrons and pions, respectively, with E in GeV. 

4 Event selection and jet definition 

A three-level trigger was used to select events online [20] . The neutral current DIS events 
were selected offline using criteria similar to those reported previously [25]. The main 
steps are outlined below. 

The scattered electron was identified using an algorithm based on a neural network [26]. 
The kinematic variables Q^, xbj and the inelasticity y were reconstructed using the double- 
angle method (DA) [27], where the hadronic final state was reconstructed using combina- 
tions of CTD tracks and energy clusters measured in the CAL to form energy-flow objects 
(EFOs) [28]. 

The following criteria were applied offline to select DIS events: 

• a scattered electron with energy E'^ above 10 GeV, to ensure a well reconstructed 
electron and to suppress the background from photoproduction events, in which the 
scattered electron escapes undetected in the rear beampipe; 

• 40 < 5 < 65 GeV, where 6 = " where Ei and Pz^i are the energy and 
Z-component of the momentum of each EFO and the scattered electron. This cut 
removed events with large initial-state radiation and reduced the background from 
photoproduction events; 

• He < 0.95, where i/e = I — 2#(1 ~ "^^^ 9 J and 9^ is the polar angle of the electron. 
Along with the previous requirements, this reduces the photoproduction background 
to a negligible level; 
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• |X| > 24 cm or \Y\ > 12 cm, where X and Y are the impact positions of the positron 
on the CAL, to avoid the low-acceptance region adjacent to the rear beampipe; 

• the Z coordinate of the vertex, ^vtx, determined from CTD tracks, was required to be 
in the range l^vtxl < 50 cm along the beam axis. This cut removed background events 
from non-ep interactions; 

• 0.04 < t/DA < 0.7; 

• 20 < Ql^ < 100 GeV^; 

• 0.0004 < xda < 0.005. 

After this selection, the jets were identified using the /cr cluster algorithm [29] in the 
longitudinally invariant inclusive mode [30] applied in the Breit frame [31] on the CAL 
and FPC cells, excluding those belonging to the scattered electron. The reconstructed 
jets were then boosted back to the laboratory frame. Jet-energy corrections were applied 
in order to account for the energy loss in the inactive material in front of the detector. 
The events were required to have at least one jet satisfying the following criteria in the 
laboratory frame: 

• the transverse energy of each jet was required to be E^r^^ > 5 GeV; 

• the pseudorapidity of each jet was required to be in the interval 2 < rf'^^ < 4.3; 

• x^'^^ > 0.036, which selects forward jets with large energy; 

• 0.5 < < 2, which suppresses the DGLAP-type evolution; 

• jets with 2.8 < rf^^ < 3.35 and with the azimuthal angle of the jet, (p expressed in 
radians, in the ranges < (f^''^ < 0.4, 1.0 < </>j"* < 2.2, 2.7 < (j^""^ < 3.6, 4.2 < (j^""^ < 5.3 
or 5.7 < (/^"^^ < 6.3 were rejected due to poor reconstruction caused by the large cell 
size in the FCAL; 

Using the sample described above, the triple differential cross sections, with the 0.5 < 
(E^tY/Q^ < 2 cut removed, were measured in two intervals of Q^, 20 < < 40 GeV^ 
and 40 < g2 < 100 GeV^ and for E}^' > 5 GeV. 

For the "forward jet-|-dijet" analysis, the events were required to have one forward jet, 
satisfying the same selection criteria as above, with the exception of the {E-I^^Y/Q'^ cut, 
and at least two additional jets with E'^* > 5 GeV. The two additional jets, chosen with 
the highest transverse energy were required to lie in the pseudorapidity region —1.5 < 
j^jet ^ 4_3_ three selected jets were ordered in 77-''^* as described in Section 2. 
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5 Monte Carlo simulations 



Various MC samples were generated. Samples of Lepto and Ariadne generators were 
used to simulate the detector response to jets of hadrons and to determine the hadroni- 
sation corrections needed for comparison with perturbative QCD calculations. The same 
samples were used to compare the measurements to the Lepto and Ariadne models. 
In addition, the data were compared to expectations of the Ariadne generator with its 
newly tuned proton-remnant treatment, "Ariadne tuned", and to the CASCADE gener- 
ator. 

The Lepto [32] MC program is based on first-order QCD matrix elements supplemented 
with parton showers (MEPS), which follow DGLAP evolution. 

The CDM approach is represented by the Ariadne 4.08 [33] MC program and its tuned 
variant ("Ariadne tuned"[f|) as described in a recent HI publication [18]. 

The Cascade MC program [19, 34] is based on the CCFM evolution and uses kx- 
factorisation of the cross section into an off-shell matrix element and an unintegrated 
parton (gluon) density function (uPDF). Predictions of Cascade were obtained with the 
J2003 set-1 and set-2 uPDFs [35]. The J2003 set-2 includes non-singular terms in the 
splitting function and reduces the cross sections at low xbj. 

In all MC models, the fragmentation of the final-state partons has been performed using 
the Lund [36] string model as implemented in Jetset 7.4 [37]. 

Both Lepto and Ariadne were interfaced to Heracles 4.6.1 [38] via Djangoh 1.1 [39]. 
The Heracles program simulates first-order electroweak radiative corrections. The 
CTEQ5L [40] proton parton distribution functions (PDF) were used in both cases. 

The events generated with Lepto and Ariadne were passed through the Geant 3.13- 
based [41] ZEUS detector and trigger simulation programs [20]. They were reconstructed 
and analysed using the same program chain as the data. 

6 Acceptance correction and systematic studies 

The correction factors to the data for detector-acceptance effects were obtained with the 
Ariadne and Lepto MC programs. These correction factors were calculated bin by bin 

^ The following parameters have been changed from their default values: the powers, defining the 
fraction of the proton remnant participating in the emission, PARA(10)=1.2 (default=1.0) and 
PARA(25)=1.2 (default=2.0); the square root of the mean value of the primordial p\ in the pro- 
ton remnant, PARA(27)=0.9 (dcfault=0.6). 
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as 

ATdet 

p _ ^''MC 
^Acc — , .had ' 

where iV^*^^ (-^mc) ^'^^ number of jets in bins of the detector (hadron) level distribution. 
For this approach to be valid, the distributions in the data must be well described by 
the MC simulation at the detector level, a condition which was in general satisfied by 
both Ariadne and Lepto. The average between the correction factors obtained with 
Ariadne and Lepto was taken. The values of Cacc were generally between 0.4 and 1.2 
for the inclusive forward-jet sample, and 0.5 to 1.4 for the "forward jet+dijet" sample. 

To ensure the correct MC reconstruction of the jets near the boundary between FCAL and 
FPC, jet profiles in the data and MC were compared and found to be in good agreement. 

The major sources of systematic uncertainty were as follows (the effects on the cross 
sections are shown in parentheses): 

• the largest uncertainty resulted from the model dependence of the acceptance correc- 
tions. This uncertainty was estimated using the deviations of Lepto and Ariadne 
corrections from their average (< 10% for inclusive forward jet sample, < 25 % for 
forward jet+dijet sample); 

• ± 3% shift of jet energies due to the CAL energy-scale uncertainty (< 10 % ); 

• ± 10% shift of jet energies due to the FPC energy-scale uncertainty (~ 15% for the 
last 77-''^* bin, negligible elsewhere); 

• the selection of inclusive DIS events (< 1%). The cuts on the scattered-electron energy, 
the X and Y position of the electron, 5 and Z^tx were varied. 

These systematic uncertainties, except for the energy-scale uncertainty, were added in 
quadrature separately for the positive and negative variations in each bin. The energy 
scale uncertainties, which are correlated between bins, are shown separately. The uncer- 
tainty in the luminosity of ± 2.2% is not included in the figures. 



7 NLO QCD calculations 

For inclusive forward jets, fixed-order calculations were performed with the DiSENT [15] 
code, at 0[ao?^ in the MS renormalisation and factorisation schemes. The number of 
flavours was set to five; the renormalisation (/xr) and factorisation (/ip) scales were both 
set to /iR = /ip = Q. The CTEQ6M [42] parameterisation of the proton PDFs was used. 
The theoretical uncertainty in the calculations was estimated considering the following 
three sources: 
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• yUR was varied up and down by a factor of two, which contributed up to 57% depending 
on the phase-space region; 

• /^F was also varied up and down by a factor of two and the resulting uncertainty was 
less than 5% except in the lowest Xbj bin and the most forward region; 

• the PDF uncertainty was estimated using the 40 different sets of CTEQ6 parton 
distribution functions [42]. 

For "forward jet+dijet" cross sections, the program NlojetH — \- [16] was used. This 
program calculates three-jet production in DIS at NLO 0{aal) and uses the MS scheme. 
The number of flavours, /ir and fip were set as in Disent. 

In order to compare the data to NLO calculations, corrections from the parton to the 
hadron level, Chad, were determined in each bin. The hadronisation corrections Chad 
were calculated as an average between those obtained from Lepto and from Ariadne 
and applied to the NLO calculations. The uncertainty of the hadronisation correction 
Chad was assumed to be the absolute difference in the two values. 



8 Results 

Cross sections for inclusive forward jets and for events containing a dijet system in addition 
to the forward jet, were measured in the kinematic region given by 20 < < 100 GeV^, 
0.04 < y < 0.7 and 0.0004 < xbj < 0.005. The differential jet cross sections are presented 
as functions of the variable ^ = Q^, xbj, -E^*, V^*- They were determined as 

C-A^ Cacc' 

where A^^^ta number of jets in a bin of width A^, C is the integrated luminosity, and 

^QED _ jV^°(?™/A''^™, where N^P (N^^^^) is the number of events selected at the 
hadron level in a given bin in the MC sample generated with (without) QED radiation. 

The triple differential cross sections were obtained in a similar manner, 
8.1 Inclusive forward-jet measurements 

The measured differential forward-jet cross sections as functions of , xbj, -Ej?* and 77-'*^* 
are shown in Fig. [2l where they are compared to NLO calculations. 
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The calculations predict lower cross sections than obtained from the data by as much as a 
factor two; however, they have a large theoretical uncertainty. The strong dependence of 
the calculation on /ir can be related to the fact that in this kinematic region higher-order 
terms become relevant. As a demonstration, the leading-order calculation is also shown 
in Fig. [2] for each differential cross section. It is far below the measurement, indicating 
that the contribution of 0{a'^) terms is significant. A recent publication [18], which 
used a harder renormalisation scale (the average of the dijets coming from the hard 
scattering), reported a smaller renormalisation-scale uncertainty. 

A comparison of the data with the Ariadne and Lepto MC is shown in Fig. [31 The 
predictions of the CDM obtained with "Ariadne default" are in fair agreement with the 
data with the exception of high i?^* and high 77-''^*, where Ariadne overestimates the cross 
sections. An investigation has shown that in Ariadne the proton-remnant fragments are 
generated with high px, therefore they show up at much lower r] than in other generators. 
The newly tuned Ariadne, also shown in Fig. [3l yields lower cross sections, in particular 
at high E^^^ and high r]^°^, and provides a good description of the data. 

The predictions of the Lepto MC are found to be in agreement with data in shape for all 
distributions, however the absolute normalisation is below the measurements by a factor 
of two. 

The measurement of differential forward-jet cross sections is compared to the prediction 
of the Cascade MC model in Fig. m Neither of the investigated uPDF sets gives a sat- 
isfactory agreement with the measurements in all distributions, suggesting that a further 
adjustment of the input parameters of the Cascade model is necessary. 

8.2 Triple-differential forward-jet cross section 

The triple-differential forward-jet cross sections as a function of 77-'''* are presented in two 
intervals of and three intervals of (i?;^*)^ in Fig. O Also shown in Fig. [S] are the expec- 
tations of the NLO calculations from Disent. The calculations generally underestimate 
the cross sections. The largest discrepancy between the data and the theory is seen in 
the high-Q^ range and for {E^^^Y < 100 GeV^. This region is sensitive to multigluon 
emission, which is lacking in the NLO calculations. 

In Fig. [6l the data are compared with Lepto and Ariadne. As was already observed in 
Fig. [3l the Lepto MC is always below the measurements. 

The cross sections of "Ariadne tuned" are below those of the "Ariadne default" in 
all the presented phase space. The difference between the two versions is smallest in the 
lowest- (E';^*)^ interval, where both are close to the data. In the highest- (E';^*)^ interval, 
the difference is big and "Ariadne tuned" gives a good description of the data. 
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A comparison of the data with the Cascade MC with two sets of uPDFs is shown in 
Fig. [3 

The expectations of Cascade are close to the data in the low-Q^ interval for set-1, while 
in the high-Q^ interval the set-2 gives a better description of the data. None of the sets 
can accomodate all the features of the data. 

8.3 Forward jet+dijet measurements 

The measured cross sections for events with a dijet system in addition to the forward 
jet are compared with fixed-order QCD calculations and LO parton-shower MC models 
in Figs. [SlfTUl The cross sections are presented as a function of Ar]i and Ar/2, and as a 
function of Ar]2 for two cases, namely, Ar]i < 1 and Arji > 1. 

A comparison between data and the predictions of NlojetH — \- is shown in Fig. [HI As 
already observed by the HI experiment [18], the NlojetH — \- calculations agree well with 
the data at large Ar]2, while they do not describe the data at small Ar)2, especially for 
small A?7i. The large Ari2 kinematics at low xbj favours dijets originating from photon- 
gluon fusion, with an additional gluon responsible for the forward jet. This case is well 
treated by NlojetH — h. The small Arji and Ar]2 region corresponds to the event config- 
uration in which all the three jets tend to go forward, away from the hard interaction. 
This configuration favours multigluon emission, which is not expected to be described by 
Nlojet-K-K. 

The comparison between data and the Lepto and Ariadne MCs is shown in Fig. |9l 
As before, the Lepto predictions are below the data for all differential cross sections. 
The "Ariadne default" overestimates the cross sections. This implies that energetic 
multiple jets are produced too often in the "Ariadne default" . The tuning of the Ari- 
adne parameters brings this model into very good agreement with data for all differential 
distributions. 

The comparison of the CASCADE MC to the data is presented in Fig. [TOl As before. 
Cascade does not satisfactorily reproduce the measurement. 

9 Summary 

A new measurement of the inclusive jet cross sections has been performed in an extended 
forward region, 2 < rf'^^ < 4.3, with higher statistics and smaller systematic uncertainties 
compared to previous studies. The measured differential cross sections are presented as 
functions oiQ^, xbj, -E^^* and 77-'''*. The measurements were compared to the predictions of 
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next-to-leading-order QCD calculations, which were found to be below the data, in certain 
regions by as much as a factor of two. The large contribution of next-to-leading-order 
corrections and the size of the theoretical uncertainty indicate that in this phase space 
higher-order contributions are important. The best overall description of the inclusive 
forward-jet cross sections was obtained by the newly tuned Ariadne MC model. The 
Cascade MC with J2003 set-1 and J2003 set-2 for unintegrated gluon density failed to 
satisfactorily describe the data. Therefore, these measurements can be used for further 
adjusting the input parameters of the Cascade model. 

The measurement of the cross sections of the events containing a dijet system in addition 
to the forward jet is presented as functions of pseudorapidity separation between jets 
composing the dijet, Arji, and pseudorapidity separation between forward jet and dijet 
system, A?72, for all Arji values and for Arji < 1 and Ar/i > 1. NLO calculations describe 
the data at large Ar]2 but underestimate the cross sections at small Ar]2, especially for 
small values of Ar/i, where, in the case of small xbj, the contribution of multiple gluon 
emission is expected to be large. The predictions of Lepto are significantly below the 
data. Ariadne with default parameters significantly overestimates the cross sections 
whereas the new tuning provides a good description of the data. The Cascade MC, 
as in the inclusive case, does not provide a satisfactory description of measured cross 
sections. 
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g2 bin 

I V.' V V 1 


da/dQ^ 
fnb/GpV^l 

\yju/ V.' c V 1 


^stat 


^syst 


ScAL 


Sfpc 


Cqed 


Chad 


20 - 30 
30 - 40 
40 - 50 
50 - 60 
60 - 80 
80 - 100 


23.21 

15.49 

9.76 

6.65 

3.21 

1.36 


±0.43 
±0.33 
±0.26 
±0.21 
±0.10 
±0.07 


+0.81 
-0.79 

+0.48 
-0.54 

+0.26 
-0.37 

+0.27 
-0.22 

+0.12 
-0.08 
+0.07 
-0.03 


+1.39 
-1.95 
+0.95 
-1.05 
+0.54 
-0.77 
+0.40 
-0.36 
+0.20 
-0.18 
+0.06 
-0.08 


+3.06 
-3.03 
+1.73 
-1.61 
+0.91 
-1.07 
+0.55 
-0.51 
+0.22 
-0.22 
+0.08 
-0.09 


0.97 
0.96 
0.94 
0.97 
0.97 
0.94 


0.89 
0.89 
0.89 
0.91 
0.89 
0.91 



Table 1: The differential cross section, da/dQ"^, in bins ofQ^ for inclusive forward 
jets. The statistical (Sstat), systematic (Sgyst) and jet-energy-scale uncertainties for 
CAL and FPC (5c al o-nd 5ppc ) are shown separately. The multiplicative correction 
applied to correct for QED radiative effects (Cqed) and for hadronisation effects 
(Chad) are shown in the last two columns. 



xb] bin 


da/dxBi 
(nb) 


Sstat 


^syst 


ScAL 


SpPC 


Cqed 


Chad 


0.0004 - 0.001 


275 


±6 


+19 
-19 


+16 
-20 


+34 
-35 


0.97 


0.86 


0.001 - 0.002 


200 


±4 


+10 
-9 


+14 
-14 


+23 
-20 


0.97 


0.88 


0.002 - 0.003 


125 


±3 


+2 
-3 


+6 
-9 


+11 
-12 


0.96 


0.90 


0.003 - 0.004 


89 


±2 


+4 
-3 


+5 
-6 


+9 
-9 


0.95 


0.92 


0.004 - 0.005 


65 


±2 


+1 
-1 


+4 
-5 


+6 
-6 


0.94 


0.95 



Table 2: The differential cross section, da/dx^j, in bins o/xbj for inclusive for- 
ward jets. The statistical, systematic and jet- energy- scale uncertainties are shown 
separately (see the caption of TableU\). 



E}^' bin 
(GeV) 


da/dE}^' 
(pb/GeV) 


^stat 


^syst 


^CAL 


^FPC 


Cqed 


Chad 


5 - 6.5 
6.5 - 8 
8 - 9.5 
9.5 - 11 
11 - 14 


245.1 
123.2 
42.04 
13.38 
2.21 


±3.2 
±2.4 
±1.40 
±0.73 
±0.21 


+6.9 
-7.6 
+4.5 
-4.5 
+1.22 
-1.20 
+1.17 
-0.99 
+0.16 
-0.20 


+15.0 
-18.7 
+7.7 
-8.1 
+1.9 
-2.35 
+0.63 
-1.49 
+0.05 
-0 


+32.6 
-31.3 
+10.0 
-10.2 
+1.50 
-1.69 
+0.62 
-1.35 
+0.08 
-0.21 


0.96 
0.97 
0.96 
0.96 
0.93 


0.91 
0.92 
0.89 
0.87 
0.89 



Table 3: The differential cross section, da /dE^^ , in bins of E-^rf. for inclusive for- 
ward jets. The statistical, systematic and jet-energy-scale uncertainties are shown 
separately (see the caption of TableUl)- 



XX 





(pb) 


5stat 


^syst 


ScAL 


Sfpc 


Cqed 


Chad 


2 - 2.3 
2.3 - 2.6 
2.6 - 2.9 
2.9 - 3.5 
3.5 - 4.3 


118.0 
361.4 
438.2 
331.0 
211.6 


±5.1 

±8.5 
±9.6 
±8.1 
±4.2 


+7.8 
-7.8 
+20.0 
-21.0 
+14.9 
-12.4 

+19.2 
-20.2 
+6.6 

-7.7 


+8.4 
-6.8 
+24.2 
-34.1 
+40.9 
-48.7 
+17.0 
-25.2 
+4.8 
-5.1 


+0 
-0.1 
+0 
-1.3 
+0.8 
-1.0 
+15.5 
-15.4 
+54.9 
-79.1 


0.97 

0.95 
0.96 
0.95 
0.98 


0.87 
0.90 
0.93 
0.92 
1.01 



Table 4: The differential cross section, da/drj^^^, in bins ofr]^^^ for inclusive for- 
ward jets. The statistical, systematic and jet- energy- scale uncertainties are shown 
separately (see the caption of TableU\). 



Arji bin 


da/dArji 
(pb) 


^stat 


^syst 


^CAL 


^FPC 


Cqed 


Chad 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.2 - 4 


82.88 
79.35 
50.68 
15.30 


±3.77 
±3.65 
±2.62 
±0.99 


+12.88 
-12.88 
+11.18 
-11.01 

+5.53 
-6.12 

+2.35 
-1.67 


+8.24 
-10.9 
+9.40 
-8.66 
+4.79 
-6.91 

+1.16 
-2.42 


+4.24 
-4.24 

+5.18 
-4.13 
+3.68 
-3.49 
+1.36 
-2.10 


0.98 
0.99 
0.98 
0.96 


0.74 
0.76 
0.76 
0.75 



Table 5: The differential cross section, da/dArji, in bins of Arji for 'forward 
jet-hdijet" events. The statistical, systematic and jet- energy- scale uncertainties are 
shown separately (see the caption of TableU\). 



Ar)2 bin 


da /d Arji 
(pb) 


^stat 


^syst 


5cAL 


^FPC 


Cqed 


Chad 


0.0 - 0.8 
0.8 - 1.6 
1.6 - 2.4 
2.4 - 3.2 


56.38 
76.23 
56.06 
19.44 


±2.64 
±3.33 
±2.82 
±1.27 


+5.26 
-4.72 
+5.23 
-6.18 
+9.33 
-8.69 
+4.46 
-3.98 


+5.42 
-5.97 
+7.44 
-9.52 
+5.18 
-7.31 
+2.12 
-2.70 


+1.35 
-1.42 
+2.51 
-3.09 
+4.44 
-4.79 
+2.69 
-2.58 


0.98 
0.96 
1.01 
0.97 


0.78 
0.79 
0.75 
0.69 



Table 6: The differential cross section, da/dAri2, in bins of Ari2 for "forward 
jet-hdijet" events. The statistical, systematic and jet-energy-scale uncertainties are 
shown separately (see the caption of TableU\). 



XXI 



A?72 bin 


da / dArji 
(pb) 


5stat 


^syst 


^CAL 


^FPC 


Cqed 


Chad 


0.0 - 0.8 
0.8 - 1.6 
1.6 - 2.4 
2.4 - 3.4 


14.86 
30.04 
29.48 
12.33 


±1.36 
±2.17 
±1.93 
±0.94 


+1.17 
-1.71 
+4.58 
-5.04 

+4.13 
-3.78 
+2.59 
-2.27 


+1.65 
-1.86 
+3.21 
-3.04 
+2.64 
-3.42 
+1.44 
-1.78 


+0 

-0.24 
+0.26 
-0.34 
+1.12 
-0.79 
+1.63 
-1.50 


1.00 
0.99 
1.02 
0.96 


0.82 
0.81 
0.78 
0.68 



Table 7: The differential cross section, da/dAr]2, in bins of Ari2 for "forward 
jet+dijet" events in the case of Arji < 1. The statistical, systematic and jet-energy- 
scale uncertainties are shown separately (see the caption of TableU^. 



Arj2 bin 


da / dArji 
(pb) 


5stat 


^syst 


^CAL 


^FPC 


Cqed 


Chad 


0.0 - 0.6 
0.6 - 1.2 
1.2 - 1.8 
1.8 - 2.8 


40.03 
46.58 
36.46 
10.27 


±2.28 
±2.59 
±2.74 
±1.03 


+5.32 
-4.81 
+1.66 
-2.71 
+7.01 
-5.80 
+2.91 
-2.87 


+3.62 
-4.12 

+3.87 
-5.97 
+4.08 
-5.49 
+0.80 
-1.44 


+1.19 
-1.14 

+1.58 
-1.76 
+4.38 
-5.49 
+1.91 
-2.51 


0.97 
0.96 
0.98 
0.99 


0.75 
0.80 
0.76 
0.71 



Table 8: The differential cross section, da/dAri2, in bins of Ari2 for "forward 
jet+dijet" events in the case of Arji > 1. The statistical, systematic and jet-energy- 
scale uncertainties are shown separately (see the caption of Table\J\). 
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5cAL 


^FPC 


Cqed 


Chad 


(GeV^) 




(nb/GeV^) 










20 < Q2 < 40 (GeV^) 




2.4-2.7 


525 ± 20_g2 


+42 
-63 




0.93 


0.88 


25-36 


2.7-3.1 


658 ± 24_24 


+52 
-90 


+ 1 
-2 


0.94 


0.91 




3.1-3.7 


463 ± 2013? 


+22 
-10 


+66 
-143 


0.98 


0.86 




3.7-4.3 


359 ± 13_i8 


+7 
-11 


+111 
-197 


0.96 


0.93 




o n o /I 


1 QQ 1 /I +10 

ioo ± 4_]^Q 


+10 
-12 




U.yo 


U.yo 




2.4-2.7 


200 ± 6+i^ 


+15 
-18 




0.98 


0.89 


36-100 


2.7-3.1 


170 ±6l^ 


+17 
-14 


+2 
-1 


0.97 


0.87 




3.1-3.7 


106 ± 4+^ 


+6 

-7 


+12 
-15 


0.97 


0.89 




3.7-4.3 


73.2 ± 2.6+1^ 


+1.7 
-1.2 


+20.9 
-27.7 


1.00 


0.88 




2.0-2.4 


17.8 ± 0.8t[!:^ 


+1.2 
-1.4 




0.99 


0.95 




2.4-2.7 


13.3 ± 0.8l[J:? 


+1.4 
-1.5 




0.99 


0.94 


100-400 


2.7-3.1 


9.56 ± QMtofj 


+0.82 
-0.95 




0.96 


0.92 




3.1-3.7 


6.14 ±0.58^°!^ 


+0.41 
-0.87 


+0.86 
-1.55 


0.97 


0.98 




3.7-4.3 


2.66 ± 0.25l[j;J^ 




+0.83 
-1.48 


0.98 


0.93 


40 < Q2 < 100 (GeV^) 




2.4-2.7 


83.6 ±5.0^1^ 


+4.6 
-7.1 


+0 
-0.2 


0.89 


0.90 


25-36 


2.7-3.1 


105 ± 5^^ 


+8 
-11 




0.94 


0.90 




3.1-3.7 


78.3 ± 4.71^2 


+2.0 
-4.5 


+7.8 
-12.0 


0.95 


0.89 




3.7-4.3 


57.2 ± 3.0+^;^ 


+0.4 
-0.9 


+15.7 

-26.3 


0.94 


0.85 




2.0-2.4 


24.2± 10.0+i;^ 


+1.5 
-1.4 




0.94 


0.92 




2.4-2.7 


35.8 ± 1.4+1^ 


+3.1 
-2.3 




0.98 


0.93 


36-100 


2.7-3.1 


26.5± 1.2ti| 


+2.2 
-2.2 


+0.1 
-0.2 


0.96 


0.88 




3.1-3.7 


19.8±0.9li;^ 


+0.9 
-0.9 


+1.4 

-1.5 


0.96 


0.89 




3.7-4.3 


13.7 ±0.6+1];^ 


+0.2 
-0.2 


+3.1 
-4.1 


0.99 


0.94 




2.0-2.4 


3.59±0.18l[J:^^ 


+0.29 
-0.30 




0.97 


0.97 




2.4-2.7 


2.57 ±0.17+°:^° 


+0.17 
-0.33 




0.96 


0.97 


100-400 


2.7-3.1 


1.71 ±0.13l[5:J^ 


+0.15 
-0.15 




0.96 


0.94 




3.1-3.7 


0.99 ± 0.08lg:J^ 




+0.13 
-0.15 


0.96 


0.96 




3.7-4.3 


0.58 ± 0.05l[5;[J^ 




+0.19 
-0.30 


0.97 


0.98 



Table 9: The differential cross sections as a function of rj^'^^ in different bins 
of and {E^'^)'^ for inclusive forward jets. The jet-energy-scale uncertainties are 
shown separately (see the caption of TableUl)- 
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Figure 1: A schematic diagram of an interaction in which a forward jet and two 
additional hard jets can be produced. 
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Figure 2: Measured differential cross sections as a function of (a) , (b) xbj, 
(c) E^* and (d) rj^'^^ for inclusive jet production (dots) compared with the NLO 
QCD calculations (solid line). The hatched area shows the theoretical uncertainties 
and the shaded area shows the uncertainty after varying the CAL and FPC energy 
scales. The inner error bars indicate the statistical uncertainties, while the outer 
ones correspond to statistical and systematic uncertainties ( except the energy-scale 
uncertainty) added in quadrature. 
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Figure 3: Measured differential cross sections as a function of (a) Q^, (b) xbj, (c) 
E^^ and (d) rj^'^^ for inclusive jet production (dots) compared with the ARIADNE 
(solid histogram), ARIADNE with new tuning (dashed histogram) and LEPTO 
(dotted histogram) predictions. The shaded area shows the uncertainty after vary- 
ing the CAL and FPC energy scales. The inner error bars indicate the statistical 
uncertainties, while the outer ones correspond to statistical and systematic uncer- 
tainties (except the energy-scale uncertainty) added in quadrature. 
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Figure 4: Measured differential cross sections as a function of (a) Q"^, (b) xbj, (c) 
E^r^^ and (d) rj^^* for inclusive jet production (dots) compared with the CASCADE 
set-1 parametrisation (solid histogram) and CASCADE set-2 (dashed histogram) 
predictions. The shaded area shows the uncertainty after varying the CAL and FPC 
energy scales. The inner error bars indicate the statistical uncertainties, while the 
outer ones correspond to statistical and systematic uncertainties ( except the energy- 
scale uncertainty) added in quadrature. 
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Figure 5: Measured differential cross sections as a function of rf^* in different 
bins of and {E^^y for inclusive jet production (dots) compared with the NLO 
QCD calculations (solid line). The hatched area shows the theoretical uncertainties 
and the shaded area shows the uncertainty after varying the CAL and FPC energy 
scales. The inner error bars indicate the statistical uncertainties, while the outer 
ones correspond to statistical and systematic uncertainties ( except the energy-scale 
uncertainty) added in quadrature. 
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Figure 6: Measured differential cross sections as a function of rj^^* in differ- 
ent bins of and {E^*^y for inclusive jet production (dots) compared with the 
ARIADNE (solid histogram), ARIADNE with new tuning (dashed histogram) and 
LEPTO (dotted histogram) predictions. The shaded area shows the uncertainty af- 
ter varying the CAL and FPC energy scales. The inner error bars indicate the 
statistical uncertainties, while the outer ones correspond to statistical and system- 
atic uncertainties (except the energy-scale uncertainty) added in quadrature. 
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Figure 7: Measured differential cross sections as a function ofrj^^^ in different bins 
of and {E^^y for inclusive jet production (dots) compared with the CASCADE 
set-1 parametrisation (solid histogram) and CASCADE set-2 (dashed histogram) 
predictions. The shaded area shows the uncertainty after varying the CAL and FPC 
energy scales. The inner error bars indicate the statistical uncertainties, while the 
outer ones correspond to statistical and systematic uncertainties ( except the energy- 
scale uncertainty) added in quadrature. 
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Figure 8: Differential cross sections for forward +dijet sample as a function of 
(a) Aril, (b) A?72, (c) A?72 for A?7i < 1 and (d) Ari2 for Arji > 1. The data (dots) 
are compared with the NLO QCD calculations (solid line). The hatched area shows 
the theoretical uncertainties. The shaded area shows the uncertainty after varying 
the CAL and FPC energy scales. The inner error bars indicate the statistical uncer- 
tainties, while the outer ones correspond to statistical and systematic uncertainties 
(except the energy-scale uncertainty) added in quadrature. 



XXXI 



ZEUS 



150 



50 







ARIADNE 

- - - ARIADNE (tuned) 
LEPTO 




(a) 






^ r-^- 




; 






• 



^150 
^ 100 



50 







1 1 1 1 1 
Energy Scale Uncertainty 
• ZEUS 82 pb ' 


1 




(b) 


♦ 
















r.r.r.r*r.r,r.r 



Aril 



Ari, 




Figure 9: The differential cross sections for forward+dijet sample as a function of 
(a) Ar/i, (b) Ari2, (c) A?72, for Ar/i < 1 and (d) Ar]2 for Ar/i > 1. The data (dots) 
are compared with the ARIADNE (solid histogram), ARIADNE with new tuning 
(dashed histogram) and LEPTO (dotted histogram) predictions. The shaded area 
shows the uncertainty after varying the CAL and FPC energy scales. The inner 
error bars indicate the statistical uncertainties, while the outer ones correspond to 
statistical and systematic uncertainties (except the energy-scale uncertainty) added 
in quadrature. 
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Figure 10: The differential cross sections for forward +dijet sample as a function 
of (a) Arji, (h) Ari2, (c) Ari2, for Arji < 1 and (d) Arj2 for Arji > 1. The data 
(dots) are compared with the CASCADE set-1 (solid histogram) and CASCADE 
set-2 (dashed histograms) predictions. The shaded area shows the uncertainty after 
varying the CAL and FPC energy scale. The inner error bars indicate the sta- 
tistical uncertainties, while the outer ones correspond to statistical and systematic 
uncertainties (except the energy-scale uncertainty) added in quadrature. 
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